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Abstract 



> 

The renormalization-group improved finite order expansions of the QCD observables have an 



unphysical singularity in the Euclidean region, due to the Landau pole of the running coupling. 

O 

Recently it was claimed that, by using a modified Borel representation, the leading one-chain term 
in a skeleton expansion of the Euclidean QCD observables is finite and continuous across the Landau 
pole, and then exhibits an infrared freezing behaviour, vanishing at Q 2 = 0. In the present paper 
we show, using for illustration the Adler-D function, that the above Borel prescription violates the 
causality properties expressed by energy-plane analyticity: the function T>(Q 2 ) thus defined is the 
boundary value of a piecewise analytic function in the complex plane, instead of being a standard 
analytic function. So, the price to be paid for the infrared freezing of Euclidean QCD observables 
is the loss of a fundamental property of local quantum field theory. 
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I. INTRODUCTION 



The renormalization-group improved expansion of the QCD physical amplitudes and 
Green functions is plagued at finite orders by the unphysical singularity of the running 
coupling (the Landau pole). For instance, the expansion of the Euclidean Adler function in 
massless QCD 

N 

V { S{Q 2 ) = T,dna n+ \Q 2 ) (1) 

n=0 

has a pole on the spacelike axis at Q 2 = A 2 , present in the one-loop coupling 

^ = ^ = AMSW < 2 > 

Here (3q = (33 — 2ri/)/12 is the leading beta-function coefficient in QCD with N c = 3 and 
rif active flavours, and Q 2 is the energy variable defined such as Q 2 > on the spacelike 
(Euclidean) axis, and Q 2 < on the timelike (Minkowskian) one. 

At finite orders, a modification of perturbative QCD leading to regular Euclidean ob- 

n n 

servables has been made in the so-called "analytic perturbation theory" pj, [2J. The basic 
ingredient of this approach is the Kallen-Lehmann representation of the Green functions in 
terms of their spectral densities. This restores the required analyticity at finite orders of the 
perturbation expansion: the powers a n (Q 2 ) of the Euclidean coupling are replaced by a set 
of functions A n {Q 2 ) which have no unphysical singularities at Q 2 > 0. 

Beyond finite orders, the observables can be defined by a summation of the Borel type. 
The Borel transform B(u) of the Adler function has singularities on the real axis of the 
w-plane |3|: the ultraviolet (UV) renormalons along the range u < —1, and the infrared (IR) 
renormalons along u > 2 (we adopt the definition of Borel transform used in [4j ) . While the 
Borel transform is, for a wide class of functions, uniquely determined once all the perturba- 
tion expansion coefficients are explicitly given, the determination of the function having a 
given perturbative (asymptotic) expansion is, on the other hand, infinitely ambiguous; not 
only due to the singularities, but because the contour of the Borel-type integral can be varied 
rather freely, without affecting the values of the expansion coefficients of the perturbation 
series. 

In Ref. joj the authors use, in two different parts of the Euclidean region, two different 
summation prescriptions of one single perturbation series. They do so by choosing two 
different integration contours of the Borel-Laplace integral, depending on the sign of the 
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running coupling. For positive coupling, a(Q 2 ) > 0, they choose the contour along the 
positive (IR renormalon) axis, 

oo 

V PT (Q 2 ) = — / e -^°°« 2 » B(it) du, a(Q 2 ) > 0, (3) 
while for negative coupling the integral is taken instead along the negative (UV renormalon) 



axis: 



— oo 

V PT {Q 2 ) = f- / e -«/(/3oa(Q 2 )) dUj a (g2) < Q ^ 

PO ^ 

This dual definition of one physical quantity has been introduced in js| because the integral 
(JSJ) is divergent for a(Q 2 ) < and vice versa, the integral (J3| is divergent for a{Q 2 ) > 0. In 
the present paper we shall examine some consequences of this definition. 

As argued in j^l], this result can be expressed also in terms of the characteristic function 
ujv(t) defined by Neubert, by using the method developed in Ref. jj]: 

oo 

V PT (Q 2 ) = Jdru v (r)a(rQ 2 ). (5) 
o 

Of course, neither of the above integrals is well-defined: in ([3]) and (jl]) the integration path 
runs along the infrared and ultraviolet renormalons, respectively. In (jSj), the integrand 
exhibits the singularity in the coupling a(rQ 2 ) at r = A 2 /Q 2 . Regulating with the Principal 
Value and taking into account the continuity of the characteristic function wp(t) at r = 1 
4j, the authors of [5| conclude from (j3J) that the contribution of the leading chain of the 
skeleton expansion of the Adler function (or other similar Euclidean observables) is finite 
along the whole spacelike axis Q 2 > and approaches a zero limit at Q 2 = 0. So, the 
unphysical Landau pole present in finite order expansions disappeared from the all-order 
summation of a certain class of Feynman diagrams. 

Of course, since the perturbative series of QCD is ambiguous, it is possible in principle 
to implement a desired property by a suitable summation prescription, or by choosing a 
different integration contour. 

The question is however whether the change of the prescription does not violate other 
undamental requirements. It turns out that this is the case for the prescription chosen in 
{J. The deficiency regards the analytic properties in the complex plane: in [5] the authors 
calculate the leading term of the skeleton expansion of the Adler function T>(Q 2 ) only on 
the Euclidean axis. In the present paper we consider the properties of this function in the 



complex energy plane. We use the fact that the characteristic function considered in [5[] is 
also the inverse Mellin transform of the Borel function B J4]. Then, using the techniques 
developed in Refs. (], [7], we calculate the Adler function in the complex energy plane, for 
the Borel prescription adopted in [5j. The aim is to see whether the prescription satisfies 
the analyticity properties in the complex energy plane required by causality in field theories 



In the next section we review a few facts about the characteristic function and the Mellin 
transform of the Borel function. In Section III we calculate the Adler function in the complex 
energy plane, using the definition for the Borel summation adopted in and in the final 
Section we discuss the analytic properties of this function. 



II. CHARACTERISTIC FUNCTION AND INVERSE MELLIN TRANSFORM 



The characteristic function ujd appearing in Eq. ([5]) was introduced in [4j , where it was 



denoted byjwo (the same notation was used in 
with Ref. 



7]). In order to facilitate the comparison 



5| we adopt the notation u)t>(t) used in this work. As shown in 
ujv is the inverse Mellin transform |l0| of the Borel function B: 

UQ+ioo 

2-ni 



, the function 



W p( T ) = _ / duB^r 11 - 1 . (6) 



uo— zoo 

The inverse relation 

B(u)= Jdruj v (r)r- U , (7) 


defines the function B(u) in a strip parallel to the imaginary axis with — 1 < Re w < 2, where 
it is assumed to be analytic. The function lot>{t) was calculated in |4( and was re-derived 
recently in Ref. |5| in the large-/?o approximation (which gives the contribution of the leading 
one-chain term in the skeleton expansion). Using ([6]), the calculation is based on residues 
theorem: for r < 1 the integration contour is closed on the right half-w-plane, and the result 
is the sum over the residues of the infrared renormalons; for r > 1 the integration contour 
is closed on the left half-w-plane, and the result contains the residues of the ultraviolet 
renormalons. The residues of the IR and UV renormalons satisfy some symmetry properties 
[5I, but their contributions are not equal. Therefore wd(t) has different analytic expressions, 
depending on whether r is less or greater than 1. Following Ref. we denote the two 
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branches of uj v by uj^ and u;^ y , respectively (in Refs. Q.BQ , oj!^ was denoted by 
and by w^). According to the above discussion, the deformation of the contour in (J6]) 
gives 

1 



u i 1 ?{t) = — r duB^u)^- 1 - duB^r 11 - 1 , (8) 

C+ C- 

where C± are two parallel lines going from to +oo slightly above and below the real positive 
axis, and 



1 

2rd 



duB(u) r" -1 - 

c, c 



duB{u) t 



u-1 



(9) 



where C'± are two parallel lines going from to — oo slightly above and below the real negative 



axis. 



14, 



The explicit expressions of uo^ and ui^ v in the large- (3q approximation [4( are 
'7 



4 



In r 



l + r) 

1 



L 2 (-t) +lnrln(l + r) 



(10) 



-U + lnr+ - + -lnr - + (l + r 



1 



r 



^ 2 (-r _i ) -lnrln(l + r 



1 



where L 2 (x) = — J X — ln(l — t) is the Euler dilogarithm. As shown in |4j, the function ^©(r), 
defined by its two branches given in fflOl) . is continuous at r = 1 (as are also its first three 
derivatives). 

In the above relations the variable r was real. However, from Eqs. ( flOl) it is clear that 
both functions uj^(t) and uj^ v (t) are in fact analytic in the r-complex plane cut along the 
real negative axis. This fact will be useful below. 



III. ADLER FUNCTION IN THE COMPLEX PLANE 

In this section we calculate the Adler function for complex values of the energy, adopting 
the choice of the Borel-Laplace integral made in [5]. Following Ref. [5| we work in the V- 
scheme, where all the exponential dependence in the Borel-Laplace integrals ([3]) and (jl]) is 
absorbed in the running coupling, and denote by A v the corresponding QCD scale parameter. 
Also, to facilitate the comparison with & we denote the complex energy variable by Q 2 
(connected by Q 2 = —s to the s-channel energy variable used in p, 12], for which s > is 
the timelike axis). 
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Let us consider Q 2 complex, first such that \Q 2 \ > Ay. Since in this case Rea(Q 2 ) > we 
use the choice (j3J) of the Borel-Laplace integral. The integral is not well-defined, due to the 
IR renormalons along u > 2//3 . As usual, we adopt the principal value (PV) prescription, 
taking 

V PT {Q 2 )= l -[V ( ~ + \Q 2 )+V^-\Q 2 )} 1 (11) 
the quantities £> (±) (Q 2 ) being defined as 



V^\Q 2 



1 



-u/WoaiQ*)) B(u) du 



(12) 



where C± are two parallel lines slightly above and below the real positive axis (these lines 
were introduced already in Eq.fjHJ)). 

Our aim is to express the integrals (fT2l) in terms of the characteristic function ujx> defined 
in the previous section. To reach this, we shall use the method developed in Ref. f| some 
time ago and used in a similar context. To make our exposition selfcontained, we shall 
expound the method here in some detail. 

As a first step, we have to pass from the integrals along the contours C± to integrals 
along a line parallel to the imaginary axis, where the representation (JTj) is valid. This can be 
achieved by rotating the integration contour from the real to the imaginary axis, provided 
the contribution of the circles at infinity is negligible. Let us first consider a point in the 
upper half of the energy plane, for which Q 2 = \Q 2 \ e l< ^ with a phase < (f> < ir. Taking 
« = Ue' 9 on a large semi-circle of radius 7Z, the relevant exponential appearing in the 
integrals (fi~2l) is 



exp 



-K 



In 



cos 9 



>sin9 



(13) 



For \Q 2 \ > A v , the exponential is negligible at large 1Z for cos# > and sin^ < 0, i.e. 
for the fourth quadrant of the complex u-plane. The integration contour defining V^~\Q 2 ) 
can be rotated to the negative imaginary w-axis, where the representation ([7D is valid. This 
leads to the double integral 



— ioo oo 

V { -\Q 2 ) = - f du [dTuo v (T)exp 



—u In 



r\Q 2 \ 



+ k 



(14) 



The order of integrations over r and u can be interchanged, since for positive <fi the integral 
over u is convergent and can be easily performed. Expressed in terms of the complex variable 
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Q , the result is 



oo 



V { -\Q 2 ) = i [dr— = /dr ^(r)a(rg 2 ). (15) 



ft J HtQ 2 /A 2 



^^ H + ?i^(W). (16) 



Consider now the evaluation of the function V^ + \Q 2 ) given by the integral along the contour 
C + above the real axis. Naively, one might think to rotate the integration contour to the 
positive imaginary axis without crossing any singularities. However, this rotation is not 
allowed, because along the corresponding quarter of a circle sin# > 0, and the exponent 
ffT3]) does not vanish at infinity for < <ft. The way out is to perform again a rotation to the 
negative imaginary u axis, for which the contribution of the circle at infinity vanishes. But in 
this rotation the contour crosses the positive real axis, and hence picks up the contributions of 
the IR renormalon singularities located along this line. This can be evaluated by comparing 
the expression (|SJ) of the function uj^ir) with the definition (|12|) of the functions D^': they 
are connected by the change of variable r = exp[— l/(P a(Q 2 ))}. It follows that T>^ can be 
expressed in terms of T>(> as 

2ni A 2 V 

AT o 2 

The relations (TTTj) . (fl5l) and (|T6|) completely specify the function T>pt{Q 2 ) for \Q 2 \ > Ay, 
in the upper half plane Im Q 2 > : 

V PT (Q 2 )= Uru v {r)a{rQ 2 ) + ^ ^ R (A 2 V /Q 2 ) . (17) 

Using the same method, the function T>pt{Q 2 ) can be calculated in the lower half of the 
energy plane, where Q 2 = \Q 2 \e 1 ^ with — n < (ft < 0. In this case, the integral along C + 
can be calculated by rotating the contour up to the positive imaginary u axis, while for the 
integration along C_ one must first pass across the real axis and then rotate towards the 
positive imaginary axis. Combining the results, we obtain the following expression for the 
Adler function for complex Q 2 with \Q 2 \ > A 2 ,: 

°° A 2 

V PT (Q 2 ) = (dru v {r)a{rQ 2 ) ±^ ^ {A 2 V / Q 2 ) , (18) 

where the ± signs correspond to ImQ 2 > and ImQ 2 < 0, respectively. We recall that the 
first term in (Tl8l) is given by the integration with respect to u, and the last term is produced 



6 



by the residues of the infrared renormalons picked up by crossing the positive axis in the 
Borel plane. 

Let us consider now \Q 2 \ < Ay, when Rea(Q 2 ) < 0. Following we use the definition 
OH) of the Borel-Laplace integral along the negative axis. In this case the integral is not 
defined due to the UV renormalons. The Principal Value prescription will be given by (1111) . 
where now the integrals are 



where C'± are two parallel lines slightly above and below the negative axis (they were intro- 
duced already in Eq. (J9])). 

We apply then the same techniques as above, by rotating the contours C'± towards the 
imaginary axis in the u plane, in order to use the representation ([7]). If the exponential (fl3l) 
decreases we can make the rotation. If not, we must first cross the real axis and perform 
the rotation. The calculations proceed exactly as before, with the difference that now one 
picks up the contribution of the UV renormalons, according to the relation (Q. This leads 
to the expression of the Adler function for \Q 2 \ < Ay 



where the signs correspond to ImQ 2 > and ImQ 2 < 0, respectively. 

IV. ANALYTICITY 

We first recall that, in confined gauge theories like QCD, causality and unitarity imply 
that the Green functions and the physical amplitudes are analytic functions in the complex 
energy variables, with singularities at the hadronic unitarity thresholds [8(. In particular, 
the Adler function D(Q 2 ) should be real analytic in the complex Q 2 plane, cut along the 
negative real axis from the threshold —4m 2 for hadron production to — oo. In perturbative 
massless QCD the first unitarity branch point lies at Q 2 = 0, but otherwise the analyticity 
in the complex plane should be preserved. This property is implemented by the well-known 
Kallen-Lehmann representation 




(19) 




oo 



(20) 




(21) 
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where TZ(s) is the perturbative part of the R e + e - ratio. In what follows we check whether the 



function T>pt derived with the Borel prescription adopted in 



5| satisfies this requirement. 



In the previous section we derived the expressions of T>pt{Q 2 ) in the upper and lower 



halves of the Q 2 plane. Since lud(t) is a real analytic function in the cut r-plane 



, i.e. 



^vij*) = ^(t)), it follows from (ITS")) and (1201) that the values in the upper and lower halves 
of the Q 2 -plane are complex conjugate to each other: T>pt({Q 2 )*) = V* PT {Q 2 ). 

We compute now the limit of these expressions when Q 2 is approaching the Euclidean axis. 
Consider first that Q 2 tends to the real positive axis from above, in the region \Q 2 \ > Ay. 
In this case we use the expression ( fl8l) . The integrand has a pole at r = Ay/Q 2 . Writing 
explicitly the real and the imaginary part of this pole, when one adds to Q 2 a small positive 
imaginary part, we obtain, for real Q 2 > A 2 : 



T> PT {Q 2 + ie) = Re 



dr ux>(r)a(TQ z 



^[K(A^ 2 )-4W)1. (22) 
Po Q 2 



We notice now that for an argument equal to Ay/Q 2 < 1, the function ujx> coincides with 
so the last term in ff22l vanishes: the imaginary part of the singularity of the integral 
is exactly compensated by the additional term appearing in ( [181 ) . 
For Q 2 < Ay, we obtain in the same way from ( |20|) 



VpAQ 2 



ie 



Re 



dr ujv{t)(i(tQ* 



" A | [(M^v/Q 2 ) - ^ V (Av/Q 2 )} ■ (23) 



A) Q 



Again the last term in this relation vanishes, since for an argument A v /Q 2 > 1 the function 
uj-d is equal to uj% v . Moreover, one can easily see that the expressions of T>(Q 2 — ie), obtained 
for Q 2 approaching the Euclidean axis from the lower half plane, differ from (122]) and fl23|) 
only by the sign in front of the last term, which again vanishes. Thus, for all Q 2 > 0, the 
functions ( |T8l) and ( l20i) approach the same expression 



V PT (Q 2 ± ie) = Re 



dr uj v {r)a{rQ 2 ) 



(24) 



This coincides with the PV regulated integral of the Cauchy type ([SD which, as shown 
in Q], is finite and satisfies the infrared freezing. Since ujv{t) is holomorphic (infinitely 
differentiable) for all r > except r = 1, the right-hand side of fl24l) has all derivatives at all 
points Q 2 > except Q 2 = A 2 , where only the first three derivatives exist 5]. This means 



S 



that ([TBI) and (|20|) define in fact analytic functions in the regions \Q 2 \ > Ay and \Q 2 \ < Ay, 
respectively. 

However, as seen from (1101) . > (r) and uu^ (r) are different analytic functions. The 
expressions (ITS]) and (1201 show that T>pt(Q 2 ) coincides with a certain analytic function in 
the region \Q 2 \ > A v , and with another analytic function in the region \Q 2 \ < Ay. So, the 
Adler function derived with the Borel representation adopted in j^] is not analytic, but only 
piecewise analytic. This is in conflict with the principle of analyticity implemented by the 
Kallen-Lehmann representation (|2T]) . 

This result implies that the infrared freezing of the Euclidean observables achieved in [s| 
has had a price. It has been possible only at the expense of analyticity, which is a guiding 
principle and fundamental property of field theory, and to which also the authors of [5[ and [91] 
repeatedly refer. The loss is not only of an academic interest: the analytical continuation is 
the only technique to obtain the Minkowskian observables from the Euclidean ones. Indeed, 
the Minkowskian quantities like the R e + e - (s) ratio are derived from the Euclidean ones by 
analytical continuation. Since the Borel-Laplace summation, as defined in does not 
satisfy analyticity, one cannot find the Minkowskian quantities in a consistent way: more 
exactly, it is impossible to define a contour, situated entirely inside the analyticity domain, 
which connects the deep Euclidean region with the low-energy Minkowskian one. 

As already pointed out, there are many different functions having the same asymptotic 
expansion in powers of the coupling. The inconsistency of the approach of refs. 5|, |9|] is not 
so much in choosing a specific definition of the integration contour, but in changing it, for 
technical reasons, during the calculation of a physical quantity. Every definition should be 
motivated physically, and so should be every variation of it. 

From the above calculations we see that the representation of the correlation functions in 
terms of the characteristic functions by using the Mellin transform technique is particularly 
suitable for examining the analytic properties of the correlation functions and performing 
their analytic continuation to \Q 2 \ < Ay and to the Minkowskian region. We mention in this 
context the calculation performed in Ref. ||]], where the function T>pt{Q 2 ) was derived for 
\Q 2 \ > A v by using the Borel prescription ([3]), and then the result was analytically continued 
to \Q 2 \ < A v . As shown in [6], the function T> PT (Q 2 ) calculated in this way exhibits a cut 
along the real axis between and Ay. Except for this unphysical singularity, the function 
is analytic in the plane cut for Q 2 > 0. Therefore, one can define in a consistent way the 
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K(s) = l/Po + limRe [° drui^M. (25) 



Minkowskian quantity 7Z by analytic continuation from the deep Euclidean region. 

We applied this method in [7], where we investigated also the question of the infrared 
freezing along the timelike axis. We used the inverse Mellin transform representation of 
T>pt(Q 2 ) derived in |6j and calculated the Minkowskian quantity 7Z(s) by analytical contin- 
uation. We find that while the fixed-order approximants TZ^ N \s) = -ImII^(s + ie) — 1 tend 
to a finite value TZ^ N \0) = l//3o in the infrared limit s — > for all iV = 1,2,3, ... (thereby 
exhibiting infrared freezing), the corresponding Borel-summed all-orders quantity lZ(s) is 
divergent with s — > 0: 

Indeed the additional integral on the right-hand side diverges like ln 2 s/s 2 for s — > 0. 
(Of course, one expects that in full QCD this divergent behavior will be compensated by a 
similar growth of similar terms in the OPE.) 

The above discussion refers only to the calculation in perturbation theory. In Ref. [s| 
the authors add to the perturbative Adler function a nonperturbative term. Although they 
use a slightly different language, based on the ambiguity cancellation in the OPE, it can be 
seen from Eq. (81) of jlj] that the nonperturbative part added to the perturbative function 
Dpt{Q 2 ) given in our relations (fT8|) and (|20|) has the form 

V np (Q 2 ) = k^M^v/Q 2 ), (26) 

where k is a real constant. Using the fact that ujx>{^\/ / Q 2 ) behaves at small Q 2 
like Q A /A.y ln(Ay/Q 2 ) [4|, one can see from the relations ( 1241) and (j26il that the sum 
DpriQ 2 ) + T^np{Q 2 ) is finite along the Euclidean axis and vanishes at Q 2 = 0. But it 
fails to be a single analytic function in the complex Q 2 -plane, being only piecewise analytic. 

In conclusion, we have shown by an explicit calculation that the Borel prescription 
adopted in [5j is in conflict with the analyticity requirements derived from causality and 
unitarity in perturbative QCD. Moreover, the simple model for the complete Adler func- 
tion proposed in [5| cannot represent the physical observable: although it is finite in the 
Euclidean region and exhibits infrared freezing, it is not consistent with the analyticity 
properties derived from field theory. 
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